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ABSTRACT:  Ab  initio  calculations  have  been  performed  to  model  the  molecular 
adsorption  of  dimethyl  methylphosphonate  (DMMP)  on  the  (001)  surface  of  the 
anatase  form  of  Ti02.  Both  the  relaxed  unreconstructed  (l  X  l)-  and  the  (4  X  l)- 
reconstructed  clean  surfaces  have  been  analyzed.  Adsorption  occurs  via  a  Ti — O— P 
dative  bond  to  a  coordinatively  unsaturated  surface  Ti  site.  In  both  cases,  one  of  the 
two  Ti— O— Ti  bridge  bonds  at  the  Ti  adsorption  site  breaks  leading  to  the  formation 
of  a  stable  Ti=0  titanyl  group.  This  species  has  not  been  reported  in  previous  studies 
of  adsorption  on  Ti02  surfaces  but  is  seen,  in  the  present  work,  as  an  intermediate  in 
the  dissociative  adsorption  of  H20. 


1.  INTRODUCTION 

One  of  the  many  applications1,2  for  titanium  dioxide  (Ti02)  is 
in  the  catalytic  or  photocatalytic  destruction,  through  oxidation, 
of  toxic  chemicals.  One  important  class  of  such  species  consists  of 
organic  reagents  containing  a  P=0  group,  which  constitute 
hazards  to  the  environment  and  to  personnel.3  Many  of  these  are 
much  too  dangerous  for  routine  experimentation;  hence,  rela¬ 
tively  nontoxic  simulants  are  used  instead,  the  most  significant  of 
which  is  (CH3)(CH30)2P=0  (dimethyl  methylphosphonate, 
DMMP). 

Due  to  the  importance  of  this  topic,  there  is  a  substantial 
amount  of  literature4  25  dealing  with  the  interaction  of  DMMP 
with  Ti02,  either  in  the  dark  or  during  photoexcitation  of  the 
Ti02  to  generate  electron— hole  pairs.  Recent  studies26,27  have 
also  investigated  the  photocatalytic  decomposition  of  actual  toxic 
compounds  on  Ti02.  On  an  OH-free  rutile  ( 1 10)  surface  in 
ultrahigh  vacuum  (UHV),  DMMP  adsorption  is  molecular  at 
room  temperature5  7  and  occurs  via  a  Ti— 0=P  dative  bond  at  a 
5-fold-coordinated  TiSc  site.  Figure  1  shows  a  model25  for  the 
corresponding  structure  on  anatase  (101)  and  also  identifies  the 
T5c,  Ti6c,  and  02c  sites.  The  Ti6c  is  6-fold  coordinated,  and  the 
02c  is  a  2-fold-coordinated  O  atom  bridging  TiSc  and  Ti6c  sites. 
For  a  given  Ti5c  site,  adjacent  02c  sites  can  be  labeled  as  “nearest” 
(N)  or  “next-nearest”  (NN)  as  illustrated  in  Figure  1.  In  addition 
to  the  dative  bond,  C — H— 02c  bonds  form  as  indicated,  based 
on  the  CH— O  distances  and  C — H— O  bond  angles,  and  have 
been  analyzed25  using  the  Atoms  In  Molecules  (AIM)  theory. 

In  contrast  to  the  UHV  experimental  data,  computational 
results  for  DMMP  on  OH-free  rutile  (110)  and  anatase  (101) 
and  (100)  surfaces24,25  find  that  dissociation  is  energetically 
favored  over  molecular  adsorption.  It  is  noted  here  that  the 
anatase  (100)  and  (010)  surfaces  are  isostructural,  as  are  the 
(101)  and  (Oil)  surfaces.  For  anatase  (101),  molecular 
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adsorption  is  exothermic  with  an  energy  of  A £ajs  =  —19.1 
kcal/mol  as  obtained  in  a  two-dimensionally  periodic  slab  (2- 
DPS)  calculation;25  whereas,  AEajs  =  —25.0  kcal/mol  is  found 
for  a  possible  first  step  in  dissociation.  Similar  results  for 
molecular  vs.  dissociative  adsorption  are  found  for  the  OH-free 
rutile  (110)  and  anatase  (100)  surfaces. 

Most  of  the  experimental  work  in  this  area  involves  infrared 
spectroscopy  and  high-surface-area  (HSA)  powders  for  which 
anatase,  rather  than  rutile,  is  usually  the  major  constituent.28 
Furthermore,  the  ( 10 1 )  face,  which  is  the  most  stable  surface  for 
a  bulk  single  crystal  of  anatase,  makes  only  a  small  contribution  to 
these  powders.  Here  the  (100)  (or  equivalently  (010))  and 
(001)  surfaces  are  dominant.  In  virtually  all  cases,  the  HSA 
powders  are  at  least  partially  hydroxylated,  and  molecular 
adsorption  on  OH-terminated  HSA  Ti02  can  be  seen  only 
below  '•“214  K  (refs  17  and  18).  Above  this  temperature,  reaction 
with  OH  occurs  easily;  although,  DMMP  does  not  react  with 
liquid  H20  under  normal  conditions.29 

The  fact  that  experiment  shows  predominantly  molecular 
adsorption  on  OH-free  rutile  ( 1 10)  but  dissociative  adsorption 
on  hydroxylated  HSA  powders,  together  with  the  computational 
result  that  dissociation  is  energetically  favored  on  OH-free 
surfaces,  suggests  that  dissociation  is  kinetically  limited  in  the 
absence  of  surface  OH.  It  has  been  proposed8,9,11,18  that  OH 
groups  promote  the  first  step  in  the  reaction  through  nucleo¬ 
philic  attack  on  the  P  atom  by  the  hydroxyl  O  atom,  leading  to 
formation  of  a  P(— O— )2  bridge  between  Ti  sites  and  either  an 
adsorbed  CH30  group  or  free  CH3OH.  Successive  steps  lead 
eventually  to  a  tenacious  PO*  residue  that  poisons  further 
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Figure  1.  Geometry-optimized  structure  (from  ref  25)  for  DMMP 
adsorbed  on  the  anatase  (lOl)  surface.  The  dashed  lines  indicate  C — 
H— O  bonds,  and  different  Ti  and  O  sites  are  labeled.  The  atom  colors 
are  blue  (Ti),  red  (O),  green  (P),  black  (C),  and  gray  (H). 

catalytic  activity.  A  recent  study30  of  the  adsorption  of 
(CH3)(0H)2P=0  (methylphosphonic  acid,  MPA)  on  OH-free 
rutile  Ti02  (110)  combining  scanning  tunneling  microscopy  in 
UHV  and  ab  initio  computation  found  that  dissociation  to  form  a 
(CH3)P(=0)(— O— Ti5c)2  bridge  and  two  OH  sites  occurs 
spontaneously  at  room  temperature.  However,  such  a  process 
may  be  less  likely  for  DMMP,  in  which  the  two  OH  groups  of 
MPA  are  replaced  by  CH30. 

The  purpose  of  the  present  work  is  to  investigate  the  adsorp¬ 
tion  of  DMMP  on  the  OH-free  anatase  (001)  surface  that, 
together  with  the  anatase  (100),  is  the  major  constituent  in  most 
HSA  Ti02  powders.  This  is  a  necessary  preliminary  to  achieving 
an  understanding  of  the  effects  of  OH  on  the  adsorption  and 
dissociation  of  DMMP  on  different  Ti02  surfaces.  The  (001) 
surface  is  also  known  to  exhibit  a  high  degree  of  photocatalytic 
activity.31  It  will  be  seen  that  molecular  adsorption  of  DMMP  on 
anatase  (001)  is  distinctly  different  from  that  on  other  Ti02 
surfaces.  On  the  relaxed  but  unreconstructed  (OOl)-(l  X  l),  as 
well  as  on  the  more  stable  (4  X  l) -reconstructed  surface, 
adsorption  of  DMMP  leads  to  an  extensive  rearrangement  of 
surface  atoms  not  seen  for  other  Ti02  surfaces. 


2.  COMPUTATIONAL  DETAILS 

All  calculations  were  based  on  2-DPS  models  and  employed 
the  CRYSTAL  06  and  09  codes,32’33  which  use  Bloch  functions 
constructed  from  localized  Gaussian  basis  sets.  Most  of  the 
relevant  details  are  given  elsewhere24,25  in  connection  with  the 
anatase  (lOl)  and  (100)  surfaces.  For  the  (001)  surface,  a  slab 
was  formed  from  a  (4  X  4)  supercell  with  a  thickness  of  4  O— 
Ti— O  trilayers,  giving  a  total  of  12  atomic  layers  and  64  Ti02 
units  per  supercell.  A  slab  of  this  thickness  is  sufficient  for 
convergence  of  the  (001)  surface  structure  and  energy34,35  and 
has  been  used  in  previous  studies36  32  of  adsorption  on  anatase 
(001).  The  need  for  a  large  supercell  will  become  apparent  when 
the  adsorption  of  DMMP  is  discussed.  DMMP  was  placed  only 
on  the  upper  surface  with  one  molecule  per  supercell,  giving  a 
coverage  of  one  DMMP  per  16  TiSc  surface  sites.  The  initial 


DMMP  conformation  was  the  one  identified40,41  as  the  lowest  in 
energy  for  the  gas-phase  molecule. 

Two  different  computational  approaches  were  employed  in 
order  to  verify  the  lack  of  a  significant  influence  on  the  final 
results.  In  the  first,  relaxation  was  done  at  the  restricted  Hartree— 
Fock  (RHF)  level,  followed  by  a  total-energy  calculation  using 
density  functional  theory  (DFT).  In  the  second,  the  entire 
calculation  was  performed  using  DFT.  In  both  cases,  DFT 
employed  the  B3LYP  functional,  the  applicability  of  which  to 
studies  of  the  adsorption  of  organic  species  on  Ti02  has  been 
discussed  elsewhere.42  During  relaxation,  the  atoms  in  the  lower¬ 
most  O— Ti— O  trilayer  were  frozen  in  the  optimized  bulk-lattice 
positions  obtained  previously25  while  all  other  Ti02  and  DMMP 
atoms  were  unconstrained.  Geometry  optimizations  used  an 
effective-core  pseudopotential  for  Ti  and  a  modified  6-31G(d) 
basis  set  for  O,  as  previously  described.24,25  The  DFT  energy 
calculations  for  the  relaxed  structures  used  the  all-electron  basis 
sets  for  Ti  and  lattice  O,  also  as  previously  described.24,25  The 
6-31G(d,p)  basis  sets  appropriate  to  the  free  molecule  were  used 
for  DMMP  since  larger  molecular  basis  sets  led  to  linear- 
dependency  problems  in  2-DPS  calculations.24,25  Due  to  the 
large  supercell  size,  all  calculations  used  a  (l  X  l)  Monkhorst- 
Pack  grid.  Other  computational  parameters  including  the  con¬ 
vergence  criteria  for  geometry  optimization,  the  definition  of 
adsorption  energy  (AEads)  and  the  correction  for  basis  set 
superposition  error  (BSSE)  are  described  elsewhere.24,25  The 
basic  approach  outlined  above  was  tested  in  several  ways  as 
described  in  the  following  two  sections. 

3.  RESULTS 

3.1.  Adsorption  on  the  (1  x  1)  Surface.  The  ideally 
terminated  (001)  surface  (Figure  2a)  is  nonpolar  and  stoichio¬ 
metric,  with  one  02c  and  one  03c  for  every  TiSc.  It  is  also 
autocompensated  since  the  excess  electron  density  on  the  TiSc  is 
transferred  to  the  electron-deficient  02c  to  give  a  closed-shell 
configuration  on  both.  However,  recent  results43  show  that 
relaxation  is  not  complete  on  the  (OOl)-(l  X  l)  surface  and 
that  there  are  occupied  surface  states  in  the  band  gap  very  close  to 
the  valence  band  maximum  (VBM).  An  alternative  termination, 
in  which  the  02c  layer  is  removed  to  give  a  surface  consisting  of 
Ti3c  and  03c  sites,  is  neither  stoichiometric  nor  autocompensated 
and  is  expected  to  be  much  less  stable. 

Table  1  shows  the  structural  changes  resulting  from  relaxation. 
These  consist  largely  of  displacements  of  O  atoms  in  the  first  and 
second  O— Ti— O  trilayers  in  opposite  directions  along  the  [100] 
axis,  together  with  an  outward  displacement  of  02c  atoms  (“Ol” 
in  Table  l)  along  the  surface  normal.35,43,44  The  TiSc  atoms 
(“Til”  in  Table  l)  also  undergo  displacements  opposite  to  those 
of  O2o  resulting  in  short  and  long  TiSc—  02c  bonds.  As  is  seen  in 
Table  1,  displacements  beyond  the  first  Ti  and  third  O  layers  are 
very  small;  however,  the  displacements  in  the  outermost  layers 
lead  to  a  large  tensile  stress,  which  is  relieved  by  the  (4xl) 
reconstruction,45  described  in  Sec.  3.3,  or  by  adsorption.36,46 

Figure  2b  shows  the  starting  structure  for  molecular  adsorp¬ 
tion  of  DMMP  on  the  relaxed  (l  X  l)  surface.  It  has  previously 
been  shown24,25  that  Ti —  0=P  bonding  at  a  coordinatively 
unsaturated  Ti  site  is  the  most  favorable  mechanism  for  adsorp¬ 
tion  and  that  bonding  to  a  methoxy  O  atom  is  negligible. 
Previous  work25  also  shows  the  formation  of  C— H— 02c  hydro¬ 
gen  bonds  and  analyzes  these  using  the  AIM  theory.  Examination 
of  Figure  2b  shows  that  initial  structures  differing  by  a  rotation  of 
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(b)  (c) 


Figure  2.  (a)  The  ideally  terminated  anatase  (001)  surface.  The  outermost  Ti  (O)  atoms  are  5-fold  (2-fold)  coordinated,  and  the  Ti  and  O  sites  referred 
to  in  Table  1  are  labeled,  (b)  The  relaxed  (00l)  surface  with  adsorbed  molecular  DMMP  before  further  relaxation.  “Short”  and  “long”  Ti5c — 02c  bonds 
(cf.  Table  l)  are  labeled,  (c)  Same  as  (b)  but  after  relaxation  following  adsorption  of  DMMP.  The  Ti— O  site  is  labeled,  and  the  dashed  lines  indicate 
possible  H-bonds  identified  on  the  basis  of  the  CH — O  distances  (2.22  to  2.55  A)  and  C — H— O  bond  angles  (139°  to  162°). 


Table  1.  Comparison  of  Calculated  Atomic  Displacements  Due  to  Relaxation  of  the  Anatase  (OOl)-(l  X  l)  Surface11 


atom 

this  workb 

this  workc 

ref  35d 

ref  44e 

Ol 

-0.14;  +0.06 

-0.17;  +0.06 

-0.19;  +0.08 

-0.09;  +0.02 

02 

-0.10;  -0.03 

1 

p 

o\ 

1 

o 

o 

<N 

O 

o 

1 

o 

1 

-0.07;  0.00 

03 

+0.07;  -0.03 

+0.10;  -0.04 

+0.11;  +0.01 

+0.05;  +0.01 

04 

+0.01;  -0.02 

+0.01;  0.00 

+0.02;  0.00 

-0.01;  +0.01 

Til 

+0.10;  -0.08 

+0.07;  -0.09 

+0.04;  -0.06 

+0.05;  -0.04 

Ti2 

-0.01;  -0.02 

o 

o 

1 

rO 

O 

O 

1 

-0.01;  +0.01 

0.00;  +0.01 

Ti3 

0.00;  -0.03 

0.00;  -0.02 

+0.02;  0.00 

Til  — Ol  short/ 

1.745 

1.738 

1.76 

Til— Ol  lon;( 

2.195 

2.207 

2.20 

Z(Til-Ol-Til)* 

147.2° 

146.7° 

146° 

All  distances  are  in 

Angstroms.  The  first  (second)  number  gives  the  displacement  along  the  [100]  ( [001] )  direction.  All  displacements  along  [010]  are 

zero.  See  Figure  2  for  definitions  of  atoms  and  axes.  Ol  is  a  2-fold-coordinated  02c  atom.  All  other  O  atoms  are  3-fold-coordinated  03c  sites.  Similarly, 
Til  is  a  Ti5c,  while  all  other  Ti  atoms  are  Ti6c.  h  Optimization  using  RHF  with  ECPs. L  Optimization  using  DFT  with  the  B3LYP  functional  and  ECPs. 
''  I) FT  plane-wave/pseudopotential  calculation;  6  O— Ti— O  trilayers;  PBE  functional.  ‘  DFT  calculation;  Bloch  functions  constructed  from  localized 
Gaussian  basis  sets;  10  O— Ti— O  trilayers;  B3LYP  functional.  These  values  refer  to  the  “short”  and  “long”  Ti5c— 02c  bond  lengths.  The  ideally 
terminated  value  is  1.936  Ain  the  present  work.  ®  This  is  the  bond  angle  at  the  surface  02c  site.  The  ideally  terminated  value  is  155.2°  in  the  present  work. 


~90c  about  the  P=0  bond  are  inequivalent  due  to  the  presence 
of  short  and  long  Ti5c-02C  bonds.  The  various  structures  differ  as 
to  whether  C — H— 02c  bonds  form  to  short-  or  long-bonded 
02c  atoms.  It  has  also  been  shown24,25  that  C — H— 02c  bonding 
tends  to  “lock”  the  DMMP  in  whatever  rotational  orientation  is 
selected  initially.  Thus  the  relaxed  structure  may  represent  only  a 
local,  and  not  a  global,  energy  minimum.  In  the  present  work, 
initial  configurations  in  which  the  DMMP  was  rotated  by  90  or 
180°  about  the  P=0  bond,  relative  to  that  shown  in  Figure  2b, 
were  considered;  although,  only  the  structure  shown  was  relaxed 


all  the  way  to  convergence.  This  geometry  appears  to  be  the  most 
stable,  but  the  differences  among  the  relaxed  energies  are  very 
small  (<2  keal/mol).  More  importantly,  the  unusual  relaxed 
structure  described  below  was  seen  for  all  three  starting  config¬ 
urations.  The  initial  structure  in  Figure  2b  shows  the  Ti — OP 
bond  at  an  angle  with  respect  to  the  surface  normal.  An  initial 
structure  with  the  Ti— OP  bond  along  the  surface  normal  relaxed 
into  the  configuration  shown  at  an  early  stage  in  the  optimization. 

Relaxation  (Figure  2c)  leads  to  a  breaking  of  the  long  Ti5c — 
02c  bond  at  the  adsorption  site,  with  the  terminal  O  atom  then 
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forming  a  Ti=0  bond  with  its  remaining  neighbor.  This  inter¬ 
pretation  is  based  on  the  Ti=0  distance  of  1.634  A,  which  is 
much  less  than  the  computed  bulk  Ti — O  distance  of  1.967  A 
(1.936  A)  parallel  (perpendicular)  to  the  [001]  axis.  The  Ti 
atoms  at  the  Ti=0  and  DMMP  sites  relax  away  from  each  other 
in  the  [  100]  direction  by  about  —0.66  and  +0.36  A,  respectively, 
relative  to  the  positions  on  the  bare  relaxed  surface.  The  02c 
atoms  bonded  to  these  Ti  sites  also  displace  along  [100]  by  about 
0.36  A  away  from  the  adsorption  site.  The  Ti  in  the  Ti— O  and 
the  remaining  02c  to  which  it  is  bonded  also  displace  outwardly 
along  the  surface  normal  by  about  0.25  and  0.18  A,  respectively. 
Most  of  these  displacements  can  be  seen  qualitatively  in 
Figure  2c.  Other  atoms  exhibit  smaller  displacements 
(generally  <0.05  A)  in  all  directions,  including  [010].  Due  to 
the  large  displacements  in  the  vicinity  of  the  adsorption  site  the 
(4  x  4)  supercell  may  not  be  sufficiently  large  to  permit  complete 
relaxation,  and  the  AEads  value  given  below  must  therefore  be 
considered  a  lower-limit  estimate  of  the  low-coverage  adsorption 
energy.  Similar  considerations  affect  AEads  for  the  dissociative 
adsorption  of  H20  (see  Section  3.2). 

The  formation  of  a  titanyl  Ti=0  site  is  an  unusual  and 
unanticipated  result.  Previously  such  species  (with  a  bond  length 
of  1.62  A)  were  proposed47  to  form  on  the  rutile  (01l)-(2  X  l) 
surface.  However,  recent  work48  has  found  that  other  structures, 
not  involving  Ti=0,  are  more  stable  for  that  surface.  The  (0001) 
surfaces  of  oxides  with  the  M203  corundum  structure  (M  =  V,  Cr, 
Fe)  have  been  shown  both  experimentally  and  theoretically  (ref  49 
and  works  cited)  to  form  M  =  O  sites.  However,  to  the  best  of  our 
knowledge  there  are  no  prior  reports  of  Ti=0  formation  in 
response  to  adsorption  on  a  Ti02  surface.  Many  other  computa¬ 
tional  studies3  ’  '  ’  of  reactive  species  adsorbing  on  anatase 
(001)  have  also  observed  the  breaking  of  a  Ti5c — 02c  bond. 
However,  the  reagents  involved  (H20,  CH3OH,  etc.)  are  all  able 
to  release  an  H  atom,  which  reacts  with  the  singly  coordinated  O  to 
form  a  Ti— OH  site  instead  of  a  Ti=0.  The  subject  of  H20 
adsorption  is  discussed  further  in  Section  3.2.  The  Ti— 0=P  bond 
is  characterized  by  Ti — O  and  O— P  lengths  of  1.962  and  1.500  A 
respectively  and  an  angle  of  163°.  The  P=0  distance  is  0.041  A 
longer  than  that  computed  for  the  free  molecule.  After  BSSE 
correction  a  DFT  adsorption  energy  of  AEads  =  —42.5  kcal/mol 
was  found.  The  correction,  AE(BSSE)  =  6.6  kcal/mol,  was 
obtained  from  previous24,25  2-DPS  calculations  for  molecular 
adsorption  of  DMMP  on  anatase  (101)  and  (100)  surfaces  using 
the  same  methods  and  basis  sets  as  those  employed  here. 

Another  issue,  which  relates  to  the  experimental  detection  of 
Ti=0  (see  Section  3.4  below),  is  whether  it  persists  after 
removal  of  adsorbed  DMMP.  It  was  found  that,  with  DMMP 
removed,  the  long  Ti5c — O  bond  does  not  reform  but  instead  the 
Ti=0  remains  intact  with  some  additional  relaxation,  consisting 
of  a  small  (<0.2  A)  displacement  inward  along  the  surface  normal 
and  toward  the  vacant  adsorption  site.  The  resulting  energy  is 
5.0  kcal/ mol  lower  than  that  of  the  relaxed  ( 1  X  1 )  surface,  which 
represents  the  exothermic  formation  energy  of  the  Ti=0. 
Subtracting  this  from  AEads  gives  an  estimate  of  —37.5  kcal/ 
mol  for  the  Ti— 0=P  bond  energy,  which  can  be  compared  with 
values  of  —21.0  and  —19.1  kcal/mol  respectively  for  molecular 
adsorption  on  the  (100)  and  (101)  anatase  surfaces.24,25  This 
trend  parallels  that  in  surface  energy35'50  (cr  =  0.90,  0.53,  and 
0.44  J/m2  respectively  for  (001),  (100),  and  (101)),  which  indi¬ 
cates  that  the  Ti— 0=P  bond  energy  decreases  with  increasing 
surface  stability  (i.e.,  decreasing  cr).  The  rutile  (110)  surface  also 
fits  this  pattern,  with  AEads  =  —17.6  kcal/mol  and  a  =  0.31  J/m2. 


3.2.  Test  Calculations.  Several  tests  were  done  to  validate  the 
procedures  described  above  and  to  demonstrate  that  the  Ti=0  is 
not  an  artifact  of  the  calculation.  Those  involving  different  initial 
configurations  for  adsorbed  DMMP  have  been  described  in  the 
previous  section.  Geometry  optimizations  for  adsorbed  DMMP 
using  RHF  vs  DFT,  as  described  in  Section  2,  were  compared, 
and  the  differences  in  atomic  displacements  were  negligibly  small 
(~5%  or  less).  The  BSSE-corrected  AEads  for  the  relaxed 
structures  (both  obtained  in  single-point  DFT  calculations  with 
all-electron  basis  sets)  were  —42.5  vs  —40.3  kcal/mol.  A  similar 
comparison  for  the  bare  (l  X  l)  surface  has  been  discussed 
already  (cf.  Table  l).  In  another  test,  RHF  geometry  optimiza¬ 
tion  was  done  for  a  different  2-DPS  model,  one  with  eight  O— 
Ti— O  trilayers,  a  (3  X  3)  supercell  and  DMMP  placed  symme¬ 
trically  on  both  surfaces  so  as  to  maintain  the  center  of  inversion 
of  the  supercell.  All  atoms  were  unconstrained  in  the  relaxation, 
and  the  same  structure  as  that  described  above,  including  the 
Ti=0  site,  was  obtained. 

Finally,  the  adsorption  of  H20  was  considered.  Several  DFT 
plane-wave/pseudopotential  calculations38'39,46'50'51  have  found 
that,  at  low  coverage,  H20  dissociates  spontaneously  on  the 
anatase  (001)  surface.  This  is  accompanied  by  the  breaking  of  a 
Ti5c — 02c  bond  and  the  formation  of  an  H-bonded  pair  of 
TiSc — OH  sites  with  AEads  in  the  range  of —32  to  —39  kcal/ mol. 
AEads  decreases  significantly  with  increasing  coverage,46'50’51 
which  may  indicate  a  repulsive  interaction  due  to  the  large 
atomic  displacements  in  the  vicinity  of  an  adsorption  site.  This 
calculation  was  repeated,  using  the  same  model  and  methods  as 
employed  for  DMMP,  beginning  with  H20  adsorbed  molecu- 
larly  at  a  TiSc  site  on  the  relaxed  (l  X  l)  surface  (Figure  3a).  The 
H20  was  adsorbed  via  a  Ti — O  bond  with  the  2-fold  axis  normal 
to  the  surface  and  with  the  plane  of  the  molecule  lying  in  the 
(010)  plane.  The  coverage  was  one  H20  per  (4  X  4)  supercell, 
i.e.,  one  H20  per  16  surface  TiSt  sites,  but  one  per  4  TiSc  sites  on 
the  [100]  axis  along  which  most  of  the  relaxation  occurs.  Early  in 
the  relaxation  (within  6  optimization  steps)  the  long  Ti5c — 02c 
bond  broke,  and  a  Ti=0  bond  formed  with  a  minimum 
interatomic  distance  of  1.650  A  (Figure  3b).  Continued  relaxa¬ 
tion  led  to  the  transfer  of  an  H  to  form  a  Ti — OH,  resulting  in  the 
same  H-bonded  pair  of  Ti — OH  sites  as  previously  reported 
(Figure  3c).  Thus  the  Ti=0  site  is  found  to  be  an  intermediate  in 
the  dissociative  adsorption  of  H20. 

3.3.  Adsorption  on  the  (4  x  1 )  Surface.  The  atomically  clean 
anatase  (001)  surface  is  known  from  experiment  to  exhibit  a  (4  X 
l)  reconstruction.52  54  Ab  initio  modeling45  shows  that  this  is 
driven  by  the  large  reduction  in  surface  energy  that  follows  from 
the  decrease  in  the  tensile  stress  intrinsic  to  the  (l  X  1)  surface. 
Figure  4a  shows  the  ad-molecule  (ADM)  structure45  in  which 
one  of  four  02c  rows  in  the  [010]  direction  (cf.  Figure  2)  is 
replaced  with  Ti03  units  to  form  three  02c  and  one  Ti4c  site  per 
(4xl)  unit  cell.  The  increase  in  energy  due  to  the  addition  of 
coordinatively  unsaturated  surface  atoms  is  more  than  offset  by 
the  relaxation  of  the  remaining  TiSc — 02c  bonds  to  “normal” 
lengths.  As  expected,  a  Ti4c  ridge  site  is  found  experimentally55'56 
and  computationally37'46  to  be  more  reactive  than  a  TiSc  terrace 
site.  In  particular,  spontaneous  dissociation  of  H20  to  form  two 
Ti — OH  groups  has  been  reported37,55  to  occur  at  Ti4c  ridge  sites 
but  not  at  TiSc  terrace  sites.  It  is  noted  that  formation  of  the  (4  X 
l)  reconstruction  requires  annealing  at  high  temperature 
(~1050  K,  ref  52)  and  might  not  occur  at  the  lower  temperatures 
typically  used  in  studies  employing  HSA  powders.  Hence,  the 
results  described  above  for  the  (l  X  l)  surface  remain  relevant. 
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Figure  3.  Stages  in  the  dissociative  adsorption  of  H20.  (a)  The  initial 
molecularly  adsorbed  state,  (b)  After  9  optimization  steps,  showing  a 
Ti— O  group  H-bonded  to  molecular  H20.  (c)  After  an  additional  29 
steps,  showing  a  pair  of  H-bond  Ti — OH  sites.  Transfer  of  the  H  atom 
occurs  at  optimization  step  number  24. 


Following  previous  studies37,46  of  adsorption  on  the  (00l)-(4  X 
1 )  surface,  the  model  used  here  is  asymmetric,  with  reconstruc¬ 
tion  on  only  one  slab  surface  and  the  other  held  fixed  with  atoms 
in  the  ideal  bulk-lattice  positions  (cf.  Section  2). 

The  surface  bond  lengths  obtained  for  the  bare  (4xl)  surface 
agreed  well  with,  but  were  about  1%  smaller  than,  those  previously 
reported.43  A  small  fraction  of  this  difference  is  due  to  the  slightly 
smaller  bulk  [100]  lattice  constant  (a  =  3.781  vs  3.786  A)  used  in 
the  present  work  The  Ti^ — 02c  bond  lengths  of  1.965  and  1.974 A 
are  close  to  the  computed  bulk  Ti — O  bond  lengths  of  1.936  and 
1.967  A.  Figure  4b  shows  the  relaxed  structure  for  adsorbed 
DMMP,  which,  similarly  to  the  (l  X  l)  result,  exhibits  the 
breaking  of  a  T4c — 02c — Ti4c  bridge  and  the  formation  of  a 
Ti=0  site  with  a  bond  length  of  1.637  A.  Since  the  Ti=0  Ti  is 
then  5-fold  coordinated,  the  stabilization  of  this  structure  is 
thought  to  result  from  the  decrease  in  the  degree  of  coordinative 
unsaturation.  A  BSSE-corrected  AEads  of  —  58.4  kcal/mol  was 
obtained.  In  agreement  with  previous  work,37  AEads  at  a  Ti4c 
ridge  site  is  significantly  higher  than  at  a  (l  X  l)  TiSc  site.  As  in 
the  case  of  the  (l  X  l)  surface,  relaxation  after  removal  of 
adsorbed  DMMP  led  to  the  continued  presence  of  the  Ti=0 


Figure  4.  Like  Figure  2  but  for  the  (00l)-(4  X  l)  surface,  (a)  Relaxed, 
bare  (4  X  l)  surface  viewed  along  [010]  direction,  i.e.,  parallel  to  the 
ridge,  (b)  Same  surface  after  adsorption  of  DMMP  and  relaxation, 
viewed  along  the  [100]  direction  (i.e.,  perpendicular  to  the  ridge).  The 
dashed  line  in  (b)  shows  a  possible  H-bond. 


site.  However,  the  energy  in  this  case  was  only  slightly  lower,  by 
1.3  kcal/mol,  than  that  of  the  bare  relaxed  (4  X  l)  surface. 

3.4.  Experimental  Observables.  Techniques  that  may  be 
suitable  for  detecting  Ti=0  are  vibrational  spectroscopy,  ultra¬ 
violet  photoemission  spectroscopy  (UPS),  and  scanning  tunnel¬ 
ing  spectroscopy  (STS).  Although  the  v(Ti=0)  stretching¬ 
mode  frequency  has  not  been  obtained  in  the  present  work, 
experiment  shows  that  the  corresponding  modes  for  V=0, 
Cr=0,  and  Fe=0  on  the  (0001)  surfaces  of  M203  oxides  all 
occur  in  the  990—1040  cm  1  (123—129  meV)  range.49  Ob¬ 
servation  of  v(Ti=0)  in  the  presence  of  adsorbed  DMMP 
would  be  complicated  by  the  strong  DMMP  absorption  band57,58 
near  1032  cm  '.  However,  the  results  discussed  above  show  that 
the  Ti=0  persists  after  removal  of  the  DMMP.  Hence,  v- 
(Ti=0)  might  be  observable  if  DMMP  can  be  desorbed  under 
conditions  that  do  not  also  eliminate  the  Ti=0.  Experiment7 
shows  that  molecular  desorption  of  DMMP,  accompanied  by 
some  decomposition,  occurs  at  about  440  K  for  rutile  (110). 
However,  AEads  for  DMMP  on  this  surface  is  only  — 17.6  kcal/mol 
(ref  25),  which  is  significantly  lower  than  on  the  anatase  (001). 
Thus,  it  remains  to  be  seen  whether  DMMP  can  be  thermally 
desorbed  from  the  (001)  surface  without  removing  the  Ti=0. 
Another  potential  difficulty  in  vibrational  spectroscopy  is  the 
presence  of  the  anatase  transverse-optic  (TO)  phonon59  at 
838  cm  .  The  width  and  intensity  of  the  TO  phonon  peak  might 
partially  obscure  f(Ti=0)  in  an  infrared  reflection— absorption 
or  a  high-resolution  electron  energy  loss  experiment. 
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Figure  5.  Total  DOS,  for  all  atoms  in  the  2-DPS  supercell,  for  the  (l  X 
l)  surface  after  removal  of  adsorbed  DMMP  and  subsequent  relaxation. 
The  labeled  features  receive  contributions  from  the  Ti=0  bond.  The 
zero  of  energy  is  at  the  bulk  valence  band  maximum. 

In  the  model  originally  proposed47  for  the  rutile  (01l)-(2  X  l) 
surface,  Ti=0  produces  a  narrow  feature  in  the  density  of  states, 
not  seen  for  other  stoichiometric  Ti02  surfaces,  at  1  eV  below  the 
VBM.  Similar  features  are  seen  in  the  computed  (l  X  l)  density 
of  states  (DOS)  shown  in  Figure  5  and  have  been  assigned  using 
the  partial  DOS  for  individual  atoms.  These  results  were 
obtained  after  relaxing  the  surface  following  removal  of  adsorbed 
DMMP,  and  the  Monkhorst-Pack  grid  was  increased  to  (4  x  4) 
for  computation  of  the  DOS.  The  main  contribution  from  the 
Ti=0  bond  is  in  Peak  1,  at  ~0.45  eV  below  the  bulk  VBM,  which 
falls  on  top  of  structure  derived  from  other  atoms.  There  is  also 
some  Ti=0  contribution  to  Peak  2,  but  a  large  part  of  the  density 
here  is  derived  from  surface  02c  atoms.  A  similar  Peak  1,  although 
much  less  pronounced,  is  found  for  the  (4xl)  surface  after 
removal  of  DMMP  (not  shown).  Here  again  the  features  of 
interest  may  be  obscured  in  UPS  or  STS  by  structure  due  to 
adsorbed  DMMP.  However,  the  same  remarks  about  the  removal 
of  DMMP,  made  in  connection  with  v(Ti=0),  also  apply  here. 

No  empty  states  localized  on  the  Ti=0  were  found  in  the 
DOS  at  or  below  the  conduction  band  minimum  (CBM,  not 
shown).  Hence,  Ti=0  might  not  easily  be  detected  using 
electron  energy-loss  spectroscopy  (ELS)  as  a  probe  of  electronic 
excitations  associated  with  this  site.  Such  empty  states  are  seen  at 
about  1.5  to  2.5  eV  above  the  CBM,  and  ELS  transitions 
involving  the  Ti=0  might  be  obscured  by  stronger  interband 
excitations. 


4.  CONCLUSIONS 

The  molecular  adsorption  of  DMMP  on  the  OH-free  Ti02 
(OOl)-(l  X  l)  and  -(4  X  l)  surfaces  was  modeled  using  ab  initio 
calculations.  The  adsorption  is  exothermic,  with  AEads  of  —  42.5 
and  —58.4  kcal/mol,  respectively.  Relaxation  leads  to  the  break¬ 
ing  of  the  long  Ti5c — 02c  bond  in  the  case  of  the  (l  X  l)  surface, 
or  a  Ti4c — 02c  bond  in  the  case  of  the  (4  x  l),  and  the  formation 
of  a  Ti=0  bond  at  a  neighboring  Ti  site.  To  the  best  of  our 


knowledge,  Ti=0  formation  in  response  to  adsorption  has  not 
previously  been  reported  for  any  Ti02  surface.  Results  presented 
here  for  H20  adsorption  on  the  (l  X  l)  surface  suggest  that  a 
Ti=0  intermediate  is  important  in  dissociative  chemisorption, 
and  consequent  Ti— OH  formation,  for  species  with  one  or  more 
labile  H  atoms.  It  is  speculated  that  other  adsorbates  (e.g., 
acetone),  which  might  form  a  dative  bond  to  a  Ti4c  or  Ti5c 
without  decomposing  to  form  Ti— OH,  could  also  produce  a 
stable  Ti=0  site.  Also,  a  well-chosen  alternative  to  DMMP 
would  not  interfere  with  spectroscopic  detection  of  the  Ti=0, 
thus  making  desorption  unnecessary.  However,  a  “strong”  dative 
bond  may  be  required  for  Ti— O  formation  to  occur.  For 
example,  a  recent  study60  of  thiophene  adsorption  on  anatase 
(OOl)-(l  X  1)  via  a  Ti5c-S  dative  bond,  with  a  A Eads  of  only  — 
5.4  kcal/mol,  shows  no  breaking  of  the  Ti5c — 02c — TiSc  bridge 
and  hence  no  Ti=0  formation.  The  implication  is  that  a 
sufficiently  strong  adsorption  bond  to  the  Ti5c  weakens  the  long 
TiSc— 02c  bond  enough  to  overcome  the  barrier  to  breaking  this 
bond  and  forming  the  Ti=0.  The  existence  of  such  a  barrier 
would  account  for  the  fact  that  Ti=0  does  not  form  sponta¬ 
neously  during  relaxation  of  the  bare  surface. 
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